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Abstract 

We present a first-principles study to understand the phenomena of interlayer exchange cou- 
pling in Fe/Nb multilayers using the linearized- muffin-tin-orbitals method within the generalized 
gradient approximation. We find that the exchange coupling oscillates with both short and long 
periodicities, which have been examined in terms of the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
model as well as the quantum well (QW) model. We have investigated the behavior of the ex- 
change coupling by artificially varying moments of Fe atoms in ferromagnetic layers. For small 
moment of Fe, the coupling shows bilinearity in the magnetic moments implying its RKKY charac- 
ter. However, at higher moments close to the bulk Fe, the saturation of long-period oscillations is 
in accordance with the QW model. Quantum-well dispersions around the Fermi level demonstrate 
that the majority-spin bands contribute largely to the formation of quantum-well states, which 
we analyze quantitatively by making use of the phase accumulation model. Our analysis indicates 
that the quantum well model gives a better description of the oscillatory behavior of the exchange 
coupling in Fe/Nb multilayers. 
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I. INTRODUCTION 



Damped long-range oscillation of the interlayer exchange coupling (IEC) as a function of 
the spacer thickness is a well known phenomenon in magnetic multilayers [l|, [3, [J. Several 
approaches have been proposed over the years in order to explain the oscillatory behavior 
of the IEC. Among these, the two prominent models are (i) the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) model 4j and (ii) the quantum well (QW) model 5[. The RKKY interaction 
stems from the spin-polarization of the intervening conduction electrons in the spacer layer. 
Bruno and Chappert 6| have shown that the exchange coupling within the RKKY theory is 
related to the topological properties of the Fermi surface of the spacer material. Since the 
magnetic atoms in multilayers are more immersed in the spin sea of the conduction electrons 
of the spacer layer near the interface, the RKKY interaction has the strongest effect at the 
interfaces p|. However, the RKKY approach is not much effective^, 6] in describing the 
correct amplitude and phase of the coupling oscillations essentially because the interaction 
between the ferromagnetic layers and the conduction electrons is not well captured in the 
RKKY approach. In the Quantum Well model, the coupling arises due to the spin-dependent 
confinement 8[] of electrons inside the spacer medium as the size of the multilayer system 
is reduced to the nanometer range. In the quantum well picture, each layer thickness of 
magnetic as well as nonmagnetic (NM) kind in the entire multilayer stack contributes sig- 
nificantly to the coupling strength implying that the IEC is not a sheer interfacial effect 9| . 
The QW states are formed near the Fermi level when spin-polarized bands are shifted away 
from the Fermi surface because of the strong magnetization on both sides of the spacer 
medium. Such states, which can be observed directly by photoemission measurements, shift 
in energy with the spacer thickness and become closely spaced when the spacer layer ap- 
pears to be sufficiently thick. Though both the models arrive at the same period of coupling 
oscillations due to their origin in the shape of the spacer Fermi surface, they differ in deter- 
mining the coupling strength 10] . It is because the RKKY model originates from the second 
order perturbation theory whereas the QW model does not make any assumption about the 
strength of the interlayer interaction. Further, QW theory predicts additional possibilities 
like oscillating density of states, quantum well dispersions etc., especially when the spacer 
layer forms a multisheet Fermi surface (FS). Such is the case in Fe/Nb multilayers, where 
Nb spacer has three sheets of Fermi surface in the (100) plane 111 ]. 
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121 ] on the ex- 



Although many computational works exist already in the literature j^, 4, 
change coupling phenomena in systems such as Fe/Cr, Fe/Au, Co/Cu etc, where the spacer 
layers are in general NM transition metals, not much theoretical understanding has been 
gained so far for Fe/Nb multilayer system. The first experimental study on the sputtered 
Fe/Nb superlattices as carried out by Mattson et al 13| reported a weak coupling with an 
oscillation period of about 9 A at room temperature. However, some ab initio band struc- 
ture studies [14]. \\\ on Fe/Nb multilayers demonstrated an oscillation period of 4.6-6.0 A. 
Neutron reflectometry data of Rehm et al 16j for Fe/Nb multilayers have earlier suggested 
an oscillatory RKKY kind of coupling for small Nb layers. In another experimental de- 
velopment, Klose et a/[l3] have shown that hydrogen charging can modify the magnetic 
coupling in these heterostructures through the alteration of the electronic structure of the 
Nb interlayer. 

Fe/Nb multilayers seem to be interesting systems to study the interlayer exchange cou- 
pling in the sense that Fe is a strong transition-metal ferromagnet (FM) while the spacer 
layer can be a superconductor (SC) at low temperatures. Thus apart from the phenomena 
of oscillatory exchange coupling between FM layers, Fe/Nb superlattices have another inter- 
esting phenomenon known as the "proximity effect", which paves the way for new sources 
of magneto-resistance with potential applications in magneto-electronics 18|] . The presence 
of the internal magnetic field in Fe layers weakens the phenomena of superconductivity due 
to the breaking of Cooper pairs. Among various FM/SC hetero-structures, where supercon- 
ductivity gets induced in the ferromagnet by bringing it in contact with the superconductor, 
Fe/Nb multilayers have been studied experimentally 161. Il9l 1201 ] to an appreciable extent. 
One such observation by Miihge et a/(20| suggests that the thickness of the SC interface 
mainly determines the actual shape of T c versus d,F e (thickness of the FM layer) curve in 
Fe/Nb superlattices. Since Nb is the highest known T c element, its interplay with Fe con- 
tinues to remain an active field in the study of such heterostructures. 

In the present work, we elucidate the coup ling: phenomena in Fe/Nb multilayers within the 
domain of the density functional theory 2l|]. Previous calculation of Shukla and Prasad 15J 
reported the IEC for the Fe/Nb multilayer system up to 7 monolayers of Nb spacer sand- 
wiched between two Fe layers. However, it could not shed light on the detailed analysis of 
the coupling phenomena due mainly to the lack of relatively large spacer thickness, which is 
necessary for such kind of study. This work thus involves FesNb m , (m=1..16) system in order 
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to understand the coupling behavior in Fe/Nb(001) multilayers. Our calculation shows that 
QW states are indeed formed in such heterostructures and the QW model gives a better 
description of the oscillatory exchange coupling in Fe/Nb multilayers. 

The organization of the paper is as follows. In Sec. II, we briefly outline the computational 
procedure adopted in the present study. Sec. Ill, under several subsections, deals mainly 
with our results that include simultaneous discussions. We finally sum up our observations 
in Sec. IV. 



II. METHODOLOGY 



Total energy calculations have been carried out for FesNb m (m=1..16) 



multilayers in the 



framework of the linearized-muffin-tin-orbitals (LMTO) method 22|, |23|, |24| within the tight- 
binding representation. The atomic-sphere approximation (ASA) is used for the potentials 



determined self consistently in the generalized gradient approximation (GGA) 25| of the 
density functional theory 211] . Various tetragonal supercells are constructed out of Fe and 
Nb monolayers (ML) where the bcc Fe layers are stacked along the [001] growth direction 
(see Ref. 15) in ferromagnetic (FM) as well as antiferro magnetic (AFM) orientations. The 
interlayer exchange coupling (IEC ), denoted by J(m), corresponds to the energy difference 
between FM and AFM configurations per unit cell structure so that 



J(m)=Eli(m) 



Ep ot (m), 



(1) 



where m is the number of spacer layers. We calculate the total energies of all occupied states 
and minimize it between the FM and AFM configurations in the self-consistent fashion for 
each Nb thickness. The average lattice parameter with reduced lattice mismatch is taken 15J 
as 3.067 A for the present heterostructures. Linear tetrahedron method has been used for the 
Brillouin zone (BZ) integration with a maximum of 840 ^-points in the irreducible wedge of 
the surface BZ. We use the same unit cell for the FM and AFM structures to obtain reliable 
energy differences between ferromagnetic and antiferromagnetic ordering. 



4 



III. RESULTS AND DISCUSSION 



A. Interlayer coupling oscillations 

Initially, we compute the IEC as a function of the Nb spacer thickness using Eq. (1) 
for FeNb m , (m = 1..7) system. The results for J(m) are in good agreement with our earlier 



calculation 



151 ] using the FP-LAPW method. This gives us confidence in carrying out the 



present LMTO-ASA based calculations, which provide a reasonable estimate of magnetic 



moments and energy differences [26|]. The IEC for Fe3Nb m (m = 1..16) configuration is then 
computed using Eq. (1), which is shown in Fig. 1. Rapid oscillations are observed up to 
9 monolayers (ML) of Nb thickness, after which the oscillatory exchange coupling becomes 
appreciably weak. The coupling changes from ferromagnetic to antiferromagnetic configu- 
ration at about 2, 7 and 10 monolayers of Nb spacer, resulting in oscillation periods of 7.7 
and 4.6 A respectively. Note that 1 ML corresponds to the interplanar thickness of 1.5335 

A. 

To better understand the oscillation periods and the coupling phenomena in Fe/Nb mul- 
tilayers, we fit the calculated variation of J(m) with m to the following asymptotic formQ] 

2 4 

J(m) = A k sin (q k m + 4> k ) /m 2 + ^ Aisin (qirn + <pi) /m 3 , (2) 

k=l 1=3 

where A^s describe the amplitudes; qk(i)'s yield the periodicities Tk(i)(=2iT / q^Ys; and 
^>fc(i)'s the phases of the /c(/)-th mode of oscillations. Our ab initio data for J(m) in Eq. 
(1) yields a reasonably good fit to Eq. (2) with the following four periods: T\ = 4.14 ML 
(~ 6.3 A), T 2 = 5.05 ML (~ 7.7 A), T 3 = 2.86 ML(~ 4.4 A), T 4 = 20.28 ML (~ 31.1 A). 
These values fall well within the previous 

Multiple periodicities arise due to existence of the multisheet Fermi surface in the Nb spacer 
layer. The Nb thickness periodicity of 7.7 A comes close to the experimental value of 9.0 
A. This discrepancy is perhaps due to pre-asymptotic effects and the difficulty of including 
the true lattice structure in our calculations. Since in the present study, the total-energy 
calculations of the IEC are limited by 16 ML of spacer thickness, which might still be in the 
pre-asymptotic region, a slight discrepancy between the experimental result and total-energy 



calculation is expected. However, t 
with the available ab initio data[13l. 



re 4.4 and 6.3 A periods are in fairly good agreement 



14 



15 



27J].The well-fitted curve of Fig. 1 suggests that 
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the interlayer exchange coupling in Fe/Nb multilayers has significant 1/m 3 dependence in 
addition to the conventional 1 / m 2 dependence given by the RKKY theory. 

If we analyze the four periodicities that are obtained upon a reasonably good fitting of 
Eq. (2), the existence of higher harmonics is observed in the coupling function. These 
harmonics (denoted by n) add terms of the form m~^ 2+ ^sin(27inm/T + phase) with i = 0, 
1..., n > 2, and T being the fundamental period. This way the fourth harmonic (n = 5) 
results in m~ 2 sin(55m + phase) with 5 = 2n/T. Now, if we express the 4.14 ML period T\ 
as 1/Ti = 1/5 + 5/27T, the effective period[8| is obtained as T e ff = 2ir/55 = Xi/(5 — T x ) = 
4.81 ML. This value is, however, a bit lower than the fit value of 5.05 ML, which corresponds 
to T 2 . The slight discrepancy may be attributed to the uncertainties in the fit by Eq. (2). 
Expressing similarly the 2.86 ML period T 3 as 1/T 3 = 1/3 + S/2ir, the second harmonic 
(n = 3) becomes m~ 3 sin(35m + phase), which yields an effective period T' e jj = 2n/35 = 
13/(3 - T 3 ) = 20.43 ML. This value is also quite close to the fit value of T 4 (= 20.28 ML) 
with a discrepancy of about 0.7%. Both the long periods T 2 and T 4 thus turn out to be the 
"Vernier" periods of the respective short period oscillations of T\ and T3. As pointed out 
by Schilfgaarde and Harrison Q], higher harmonics have significant presence in the quantum 
well limit where the RKKY description does not hold good. 



B. Fermi surface and the RKKY periods 

In the RKKY approach, the oscillatory periods of the interlayer exchange coupling are 
uniquely determined by the stationary spanning vectors of the bulk Fermi surface of the 
spacer material. Several spanning vectors in the FS give rise to multiple periodicities, as we 
have already come across in the previous section. 

In Fe/Nb multilayers, the spacer layer Nb has five conduction electrons per atom that fill 
the first Brillouin zone completely while the second and third Brillouin zones partially 28 1. 
The second zone is a closed "octahedron" (OCT), which contains a hole sheet centered at T. 
However, the third zone has two sheets. One sheet contains an open surface of holes, referred 
to as "jungle gym" (JG), which extends from r to H points along the [100] direction. The 
other sheet is a set of distorted hole "ellipsoids" (ELL) centered at N points. The OCT 



and JG sheets contact at three points 
Fermi surface of Nb. 



291 ] in the (100) and (110) symmetry planes of the 
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Simple square lattice planes of Fe and Nb with primitive translations of -[100] and 
^[010] are stacked along the [001] growth direction to form Fe/Nb multilayers. Since the 
translational symmetry in the growth direction is broken, multilayers behave like quasi-two- 



dimensional systems, which are periodic only in two dimensions [301] . Hence, the Brillouin 
zone needs to be constructed in 2D to deal with the in-plane coordinates 3J|. Fig. 2 shows 
the cross section of the spacer layer Fermi surface of Nb on the basis of our self-consistent 
calculations in the central (100) plane. The line Qi, Q2 and Q3 as shown in Fig. 2, represent 
the spanning vectors for the [100] crystalline orientation; Qi spans the T-centered octahedron 
along the [100] direction, while Q2 and Q% span the outer and the inner ellipses respectively 
along the [100] direction. We find that the RKKY periods of 4.1, 6.4 and 7.5 A, as predicted 
from the Fermi surface spanning vectors, turn out to be in good agreement with the interlayer 
coupling periods of 4.4, 6.3 and 7.7 A, as obtained by fitting the self-consistent results to 
Eq. (2). However, the Fermi surface topology of Nb is not able to predict the periodicity 
of 31.1 A. As we know from the preceding section, the long- wavelength coupling period of 
31.1 A originates from the higher order terms in the coupling function of Eq. (2), which is 
also a higher harmonics of the short-wavelength period of 4.4 A . 

C. Magnetization 

To examine how the Fe magnetic moment at the interface behaves as a function of the 
intervening layer thickness in Fe/Nb multilayers, we plot the magnetic moment of the in- 
terface monolayer of Fe vs. the spacer layer thickness of Nb in Fig. 3. There is a good 
agreement between our calculated results and the experimental observations as shown in 
the inset. We find that the Fe magnetic moment gets saturated in the asymptotic region 
with reduction of about 25 % of the bulk value. This reduction in the measured data is 
about 40 %. Following the works of Holmstrom et al 32|, this kind of discrepancy may be 
attributed primarily to the interface alloying. 

The induced magnetic moments in the Nb spacer layer for a ferromagnetically ordered 
Fe3Nbi6 superlattices are shown in Fig. 4(a). On the other hand, Fig. 4(b) displays 
the induced polarization for the antiferromagnetic configuration. Fe layers are at positions 
marked by and 17 (not shown in the figure). The induced magnetic moment (in absolute 
values) in Nb spacer layer decreases from about 0.28 /xg at the interface to about 0.02 
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further apart from the interface. For an Nb layer of 16 atomic planes, the calculated induced 
moment at the center of the Nb layer is about 10 -3 \ib- The period of oscillation of the 
induced moment in both the magnetic configurations turns out to be about 4.6 A, which is 
in fairly good agreement with the coupling period of 4.4 A, as determined by Eq. (2). The 
bias in the magnetic ordering of the induced moment (see Fig. 4), which occurs only in the 
ferromagnetic configuration, may be regarded as being due to the onset of non-RKKY terms 



in the coupling function when the Nb layer gets appreciably thicker. Mathon et aZ[33[] have 
already shown analytically the presence of such non-RKKY terms for Co/Cu(001) system 
using the stationary phase approximation. 

However, to get a better understanding of the coupling behavior, one needs to examine 
the influence of Fe magnetization on the coupling strength. For this, we adopt a procedure 
similar to that of Schilfgaarde and Harrison^]. A trial density is constructed out of the 
charge densities of self-consistently calculated bulk ferromagnetic Fe and paramagnetic Nb in 
their respective atom-centered spheres. The charge density of Fe is constructed as follows 8|] 

/x n/ x nUr) — nHr) . . 

n Fe (r)=n°(r)±a^^ (3) 

where n° denotes the density of bulk paramagnetic Fe while n) and represent the spin 
densities of the majority and minority spins in ferromagnetic Fe and a is a parameter ranging 
from to 1. The magnetic moment of Fe atom will be proportional to a, taking full value at 
a = 1 and value when a = 0. Since we are interested in obtaining J as a function of the 
Fe moment, we have used the trial densities in our frozen-potential calculations pertaining 
to relevant supercells. To have a preliminary idea of how the coupling and the moments 
are inter- related in Fe/Nb systems, we construct a 16-atom Nb supercell with two Fe atoms 
substituted such that one Fe sits at (0 0) while the other at 3a(l 1 l)/2. For parallel and 
antiparallel alignment of two inequivalent Fe atoms, the energy difference is calculated}^ 
as a function of the Fe moment parametrized by a. According to the RKKY theory, this 
energy difference should vary as a 2 (see Fig. 5). 

Several superlattices out of Fe3Nb m (m = 1-16) multilayer configuration are con- 
structed subsequently to calculate the energy difference J a (m) = £[FejNb m Fe^Nb m ] - 
-EfFegNbmFegNbm], as function of a. We determine the amplitudes A^'s for each a by 
fitting J a (m) to the functional form as given by Eq. (2). This way, we find two short-period 
amplitudes and two long-period amplitudes, which are illustrated in Fig 5. In Fe/Cr 
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multilayers, Schilfgaarde and Harrison 8j found the coupling amplitudes (for small a) of 
both short and long periods to increase as a 2 , while at larger a, short period amplitude 
continued its rise, though the long period saturated completely. Here we find that the 
initial a 2 dependence of the exchange coupling, as assumed in the RKKY theory, at small 
a values is followed by two distinct features at higher a values viz. 

(i) stronger dependence on the Fe moment as IEC shoots up for short period oscillations 
and 

(ii) complete independence of the Fe moment as IEC gets saturated for long period 
oscillations. 



As a function of magnetization (/x), the interlayer exchange energy can be described as 35| 



J(A = J\{l*X-l**) + ■M^l-Ma) j ( 4 ) 

where J\ and J2 are respectively the bilinear and biquadratic coupling constants, while [i\ 
and fj,2 denote the magnetizations of the adjacent ferromagnetic layers. In order to fit the 
calculated results for the 16-atom supercell data, we need to expand the above expression 
in the following form 

Ja = J\ot 2 + J2« 4 + (5) 

where /j. has been parametrized by a and J3 is an additional triquadratic term appearing 
in the exchange coupling. Initially, we calculate the coupling energy, as a function of a, for 
the 16-atom Nb supercell with only two Fe atoms substituted. These are represented by 
asterisks in Fig. 5. While fitting the data we find that the coupling energy has biquadratic 
and triquadratic terms in addition to bilinear terms in a. In doing so, the coupling constants 
turn out to be in the ratio of about 1:2:4. It may be noted that the RKKY theory assumes 
only the bilinear terms in a. The linear part of the fitted curve in Fig. 5 thus shows the 
RKKY kind of coupling in the region of small magnetic moment of Fe. We then calculate 
the interlayer exchange coupling for the FesNb m (m = 1-16) multilayer configuration as a 
function of a. Short period amplitudes of the Fe3Nb m (m = 1-16) multilayer configuration 
are found to follow the above phenomenological expression for the interlayer interaction, 
but the long period amplitudes show a saturation behavior at large a. For a complete 



confinement, the coupling becomes independent of the size of the magnetic moment 36j, as 



demonstrated by the long period amplitudes. The saturation of the long-period amplitudes 
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with the ferromagnetic Fe moment thus favors a QW description of the exchange coupling. 
Because of the incomplete confinement of states in Nb, Fe/Nb multilayers exhibit a partial 
transition to the QW character from the RKKY description, which is valid at small moments 
[Fig. 5]. 



D. Energy bands, density of states and the quantum well model 

To study the nature of the electronic states in Fe/Nb multilayers, we resort to the bulk 
band structures (see Fig. 6) of the ferromagnetic Fe (for both spins) and the paramagnetic 



Nb along the [100] as well as the [110] direction. These directions are important [11] for 
multiple scattering, since the respective A2, Si and £2 bands cross the Fermi level {Ep) in 
the spacer medium contributing to the oscillations in the interlayer coupling. 

The Fe minority r 12 -A2-H 12 band closely resembles the corresponding Nb band indicating 
that the Fe/Nb interface will be more transparent to the minority-spin electrons than to the 
majority-spin electrons. The vanishing overlap between the majority spin bands Fi2-A2--£/i2 
in Nb and Fe creates spin-dependent gaps at the Fe/Nb interface resulting in energy barriers 
that confine electrons to the intervening Nb layer. The electrons in Nb band ri2-A2(T)--Hi2 
can not pass into the Fe majority band since there are no available states of A 2 (|) symmetry 
from T 12 (Fej) to T 12 (Nb). As a result, spin-up electrons with energies from (Ep - 0.14) 
eV to (Ep + 3.15) eV experience multiple reflections at the interface and get confined to 
the Nb spacer. Similar phenomena happen along the [110] direction. The electrons in the 
Nb band r 2 5'-Si(t)-A r 1 / can not pass into the Fe majority band in the energy range from 
(Ej? - 0.94) eV to (E F + 2.38) eV. Similarly, the Nb electrons in r 25 ,-E 2 tf)-N 2 band can not 
pass into the Fe majority band in the energy range from (Ep - 2.29) eV to (Ep + 0.56) eV. 
The confinement of electrons thus occurs in the spacer layer resulting in the formation of 
quantum well states. As Fig. 6 shows, bcc Fe has mostly minority-spin states at the Fermi 
level and this gives rise to majority-spin quantum well states in Nb. 

Fig. 7 displays how the density of states at Ep oscillates with the spacer layer thickness, 
which is a characteristic feature of the quantum well model. According to the QW model, 
the oscillatory behavior of the magnetic coupling in Fe/Nb multilayers as demonstrated by 
Fig. 1 can be traced back to the oscillations of the density of states (DOS) at the Fermi 

nn 

levelp, l27[. During the formation of QW states, Fe/Nb interfaces act as electron mirrors 
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that induce standing waves in the Nb spacer medium. As Fig. 7 suggests, the first two 
maxima correspond to the oscillation period of about 4.6 A. The periodicity then enhances 
to about 6.1 A. We notice that the QW periodicities of 4.6 and 6.1 A are in fairly good 
agreement with the interlayer coupling periods of 4.4 and 6.3 A as obtained from the self- 
consistent results for Fe3Nb m (m=1..16) multilayers (see Sec. III-A). It implies that there 
is a reasonable connection between the existence of the spin-polarized QW states and the 
manifestation of magnetic coupling in Fe/Nb multilayers. 

Oscillations in the density of states at the Fermi level can also be explained on the basis 
of the bulk band structures of Fe and Nb, assuming the QW model. We know from the 
previous discussion that the majority-spin energy bands of Fe along the [100] and the [110] 
directions provide the confinement to Nb electrons in forming the QW states inside the 
spacer medium. Fig. 6 shows that A2 and Si bands lead to the formation of QW states 
around Ep inside the Nb layer. It is because there are no propagating states of the same 
symmetry at the same energy associated with the majority-spin electrons upon crossing the 
Fermi level. We obtain a band offset of 4.1 eV for the A2 band in the [100] direction, and 
another band offset of 3.3 eV for the Si band in the [110] direction [see Fig. 6]. Ortega et 
al 37|] have shown for Co/Cu and Fe/Cu systems that a large band offset between the band 
edges of the magnetic and the spacer layer gives rise to sharply confined wavef unctions. 

For z being the multilayer growth direction, the translational symmetry holds only in the 



x-y plane so that the wavefunction of an electron in a quantum well can be expressed as 38] 



i/> K , k {R) = Af k (z)e lK -Pu^ k (R), (6) 

where A is a normalization constant; n = (kx,ky); p = (x,y); u Ktk (R) is periodic in the 
lattice while fk(z) is the envelope function, which ensures that the boundary conditions 
are met at the interfaces. Eq. (6) shows that the wavefunction of a QW state consists of a 
rapidly oscillating Bloch function, which is modulated by a slowly varying envelope function. 
For a single band edge, the modulation of the Bloch wave [k e dge] by an envelope (k env ) yields 



a total wavevector 



37f] 



ktot k e dg e i k envy (7) 

where only the normal components of the wavevectors are considered since the boundary 
conditions for the components that are parallel to the interface mimic the bulk. By drawing 
an analogy with the simple interferometer model, we may expect the interference maxima 
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to appear at every half wavelength (X env ) of the envelope as the thickness of the Nb spacer 
layer is increased so that QW states appear with a period of X env /2. For QW states at the 
Fermi level, the total wavevector becomes the Fermi wavevector if the spacer layer is thick 



enough to exhibit bulk-like bonding|39j. As Fig. 6(a) shows, the band edge for the Ei band 
in the [110] direction is located at the zone boundary (ksz) so that k env = [ksz - kp). The 
frequency associated with this k env in real space may be viewed as a beat frequency between 



the Fermi wavelength and the atom spacing 



101 ] . The oscillation period thus becomes 



T = Ti / k env = n I (k BZ - h F ) (8) 

This way, we obtain a period of 5.8 A due to the Ei band, which is in fairly good agreement 
with the periodicity of 6.1 A, as appeared in the oscillating DOS at Ep (see Fig. 7). It may 
be noted that the wavevector determining the periodicity of the envelope function along the 
[110] direction is quite small and thus the wavelength is large. However, a reverse situation 
occurs for the A2 band in the [100] direction, where {ksz - kp) is large as kp is small. The 
envelope function should produce in this case a short-wavelength weak modulation of the 
3- 



QW wavevector [3£ 



E. Phase accumulation model 

A simple way of predicting the energy for QW states is the phase accumulation model 
(PAM), which is often used to describe the thickness dependence of QW energies and reso- 



nances in thin films and other layered structures [401] . According to this model, the condition 
for a QW state is determined by summing over all the phases accumulated by a propagat- 
ing plane wave inside a quantum well. The total phase accumulated must be an integral 
multiple of 2n. For m monolayers of Nb spacer, each of thickness d, the total spacer width 
becomes dm = md. Thus the distance traversed by an electron in the spacer medium is 
2md resulting in 2mdk ± phase change, where k L represents the electron wavevector normal 
to the layers. If the phase shift of the electron wavefunction upon reflection at each interface 



is 0/, the total phase in a round trip within the spacer medium can be written as [40] 

2k ± d Nb + 2(j) I (E) = 2mr, (9) 
where n is an integer related to the number of half-wavelengths that span the quantum well. 
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The phase shift at each interface in Eq. (9) can be calculated by making use of following 



ansatz 



40] 



E-E L 



ME) = 2sm- 1 ^ Eu _^ -7T, (10) 

where Ejj and El represent the energies of the potential well taken from the upper and 
lower edges of the energy gap in the Fe majority-spin band, which acts as a potential energy 
barrier for the propagating electrons. The interface reflection phase evolves across a band 
gap from to -it in traversing from the top to the bottom of the potential well. 

In order to obtain k L in Eq. (9), we make use of a two-band nearly free electron (NFE) 
model, which approximates the Nb electrons in the A 2 and Si bands. According to this 
model, the wavevector of a NFE band can be expressed asf^. l41| 

k L /k BZ = l-[l + (E- V )/G - y/A(E-V )/G+(U/Gy] ^ , (11) 

where G = h 2 k BZ /2m* , where m* is the electron effective mass; Vo is a constant offset of the 
periodic potential; 2U is the energy gap at the zone boundary; and E is the electron energy 
with respect to the Fermi level. Eq. (10) thus needs three parameters to determine k L for a 
given energy E. A fit to the self-consistently calculated A 2 band of Nb yields U = 2.05 eV, 
Vo = -9.85 eV and m* = 1.08m e , where m e is the electron mass. On the other hand, upon 
fitting Eq. (11) to the Si band of Nb, we have m* = 1.05m e , U = 1.66 eV and Vq = -5.5 
eV. 

The presence of a periodic atomic potential in the spacer layer, however, ensures that in 
addition to two traveling waves with wavevectors k and -k arising from reflections at the 



interfaces, electron wavefunctions in the phase accumulation model should be described 4QJ, 
421 ] by the combination of two more waves with wavevectors (k-g) and -(k-g), where g is the 
reciprocal lattice vector. These additional waves correspond to the Bragg scattering within 
the periodic potential of the spacer layer. Since the reciprocal lattice vector normal to the 
interface within the first Brillouin zone is g = 2k BZ = 27r/d, we have 2k BZ dm = 2vmr . Using 
this relation, Eq. (9) can be rewritten as 

2K ± d Nb -2(j) I (E) = 2isir, (12) 

where fir 1 = k B z - k 1 - and v — m - n. It follows from Eq. (8) of the preceding section that 
k is characteristic of the envelope function that modulates the QW wavefunction. 
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As we know, states within Nb layers that do not coincide in energy and momentum 
with those within Fe layers form the quantum well states. With the increase in the Nb 
layer thickness, the positions of the spin polarized QW states vary and thus exhibit regular 
dispersion through the Fermi level. For n — 1, Figs. 8(a) and 8(b) demonstrate the nature 
of QW dispersions in Fe/Nb multilayers along [100] and [110] directions as generated by Eq. 
(12) of the phase accumulation model. We find that the A2 band yields weak dispersion 
as compared to the Si band of Nb due mainly to the weak modulation of k 1 - in the [100] 
direction. The oscillation periods of the quantum well states at the Fermi energy, as shown 
by the dispersion curves of 8(a) and 8(b), are obtained as about 5.5 and 7.3 A. These values 
are close to the fitted values of 6.3 and 7.7 A respectively. 



IV. CONCLUSIONS 



In this paper, we have explained the phenomena of interlayer exchange coupling in 
Fe/Nb(001) multilayers in terms of RKKY as well as the QW model using the density 
functional calculations. The RKKY periodicities arising out of extremal spanning vectors 
of the multisheet Nb Fermi surface have been found to be in favorable agreement with the 
available experimental as well as calculated data. In the region of small magnetic moments of 
Fe, both the short and long periods display the RKKY kind of coupling, since IEC shows bi- 
linearity in magnetization. However, the presence of additional biquadratic and triquadratic 
coupling constants in IEC at higher Fe moments signify the onset of the non-RKKY char- 
acter, especially of the short-period oscillations. We have found that at moments closer 
to the bulk value, the long-period oscillations become eventually independent of the size 
of the magnetic moments. The appearance of higher harmonics in the well-fitted envelope 
of the oscillatory exchange coupling turns out to be responsible for the saturation of long 
periods. The oscillatory nature of the density of states at the Fermi level supports the QW 
description of the exchange coupling in Fe/Nb multilayers. It is because such oscillations 
originate mainly from the spin-dependent confinement of the propagating electrons inside 
the quantum well of the spacer medium. Quantum well dispersions around the Fermi level 
illustrate that the majority spin bands contribute largely to the formation of QW states, 
which is subsequently analyzed by making use of the phase accumulation model. All these 
results show that magnetic quantum wells are formed in Fe/Nb multilayers and the oscil- 
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latory behavior of the exchange coupling in Fe/Nb multilayers is better described by the 
Quantum Well model. 
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Figure Captions 



FIG. 1: The calculated oscillatory interlayer exchange coupling (solid circles) as a function of the 
number of Nb spacer layers in FesNb m (m = 1-16) multilayers. The solid line is the fitted plot (see 
the text for details). The thickness of 1 Nb layer=1.5335 A. 

FIG. 2: (Color online) Cross sections of the Fermi surface of Nb in the (100) plane. T labels the 
center of the Brillouin zone, iV indicates the center of each face of the dodecahedron and H labels 
the corners of the four-fold symmetry on the zone boundary. 

FIG. 3: Variation of the Fe magnetic moment with Nb spacer thickness in Fe3Nb m (m = 1-16) 
heterostructures (see the text for details). The inset shows the experimental results by Mattson et 
alQ. 

FIG. 4: The induced magnetic moment in Nb spacer layers for Fe3Nbi6 heterostructures in the 
(a) ferromagnetic (FM) and (b) antiferromagnetic (AFM) orientations. The thickness of 1 Nb 
layer=1.5335 A. 

FIG. 5: (Color online) Interlayer exchange coupling in Fe/Nb multilayers as a function of Fe 
magnetic moment parametrized by a (see the text for details). The asterisks represent the results 
for the 16-atom supercell while the solid line depicts the fitted curve according to Eq. (5). Note 
that the linear part of the solid line resembles the RKKY kind of coupling. The triangles and 
diamonds indicate the short-period amplitudes while the circles and squares represent the long- 
period amplitudes of the FesNb m (m = 1-16) multilayer system. 

FIG. 6: The bulk energy bands of Nb and bcc Fe (| and j) along the (a) [100] and (b) [110] 
directions. The bands with A2 and Si symmetries are displayed by solid lines while the bands 
with Si symmetry by dashed lines. Only majority-spin states in Nb exhibit quantum well character 
at the Fermi level since the minority-spin A2 and Si states couple with the corresponding states 
in Fe. 

FIG. 7: Oscillations in the density of states at the Fermi level, Ep, with the Nb spacer thickness, 
caused by the quantum well states in Fe3Nb m (m = 1-16) heterostructures. 

FIG. 8: The thickness dependence of the QW energies in Fe/Nb multilayers generated by Eq. (12) 
of the phase accumulation model with respect to (a) A2 and (b) Si bands. 



